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摘要 :植物 叶片 是 碳水 交换 和 能 量 平衡 过 程 最 重要 的 场所 ,是 农林 生产 经 营 中 的 模型 佑 算 以 及 物种 结构 变异 :功能 适应 机 制 分 
析 的 关键 参考 量 。 采 用 游标 卡尺 和 手持 叶 面 积 仪 ,测量 杉木 ( Cunninghamia. lanceolata ) Župi K (LL) B KW $8 (LW,) .最 大 
WEERCLT, )3 个 直 测 指标 ,和 叶 面 积 (14) 平均 叶 宽 (ZW, ) SETWEERCLT, ) 、 叶 延长 率 ( LE) 和 叶 周 长 (LP)5 个 间接 转 算 指 
标 。 分 析 8 个 形态 学 指标 的 统计 分 布 及 其 相关 性 ,用 多 变量 线性 回归 模型 和 非 线 性 回归 指数 模型 对 7 个 形态 学 指标 和 杉木 单 
叶 叶 面积 进行 拟 合 ,结果 表明 :杉木 单 叶 面积 大 部 分 值 (95% CT) 分 布 在 0.758 一 0.836 cm” ;其 叶 面 积 的 变异 程度 最 大 (CV= 
0.513) , 叶 长 . 叶 宽 与 叶 面 积 相 关 性 达到 极 显 车 (r=0.896, 0.682), META 的 多 元 线性 模型 为 :7= -0.388+0.165 -0.023X, + 
1.453X, (R^ 20.981, SE=0.053) ,X,—X, 4194 LP LE IW garo M fi ft EBJE, LL 的 单 变量 指数 模型 适合 对 LA 进行 估算 : 
LA=0.1x(1+LL) °” (R? =0.77,X 20.39) 。 研 究 结 果 为 准确 佑 测 其 他 叶片 功能 性 状 指标 提供 了 方法 ,为 杉木 叶 面 积 估算 模型 提 
供 了 基础 数据 。 
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Abstract: Leaves are considered a significant site of carbon-water exchange and the energy-balance process. The 
characteristics of leaf response to the external environment, functional traits, coordination mechanisms, and trade-off 
strategies, as well as the structure and variation of leaves, have recently attracted huge interest. Further, leaf size could 
directly affect their capacity for light interception and carbon acquisition. Leaf area ( LA) and related leaf traits such as 
specific leaf area ( SLA ), leaf area index ( LAI), and normalized difference vegetation index ( NDVI) are the key 
indicators in crop breeding, agroforestry production and management , model estimation, and species structure variation and 
functional adaptation mechanisms. The leaf traits are based on leaf area and influenced by leaf morphology and size, and 
determination of leaf area 1s the basis to discuss plant photosynthetic production and the physiological-ecological process. 


However, the uncertainty of needle leaf area, which is due to the difficulty of measurement, could be a hindrance to 
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efficient production and management, effective risk assessment, and development of correlational research. Thus, it is 
significant to explore the appropriate measurement methods to obtain accurate values of leaf area. There is lack of 
comparative studies between different methods of measuring leaf area, resulting in inconsistency of associated concepts and 
definitions. Presently, instrument-based measurement methods of LA are pervasive and prevalent, but lack calibration by 
artificial measurements. Additionally, studies with regression estimation models based on plant leaf area and morphometric 
characteristics are concentrated on agronomy products such as crops and fruit and some broadleaf species. Research. that 
centers on needle leaf area estimation models based on leaf morphological characteristics is still lacking. This study used 
Cunninghamia. lanceolata, a common pioneer tree species in southern China, to measure 3 leaf morphometric characteristics 
(leaf length, leaf maximum width leaf maximum thickness) and indirectly calculate 5 indicators (leaf mean width, leaf 
mean thickness, leaf area, leaf elongation, leaf perimeter) using Vernier calipers and portable leaf area meter. We present 
the statistical distribution and correlation analysis of the 8 morphological characteristics, fitting the leaf area with the 7 other 
indicators in multivariate linear regression models and nonlinear regression index models. We find that (1) by manual 
measurement, the credible simple leaf area of Chinese fir ranges from 0.758 cm^ to 0.836 cm^, and shows a maximum 
coefficient of variation (CV = 0.513) ; (2) leaf area significantly correlates with leaf length and width (r- 0.896, 0.682) ; 
(3) the multivariate linear regression model of leaf length and width that is most accurate; 了 = —0.388 + 0.165X,- 0.023X, 
t 1.453X,(R* = 0.981, SE = 0.053) , where X,, X,, and X,are leaf perimeter, leaf elongation, and leaf mean width, 
respectively. From the point of view of simplicity, leaf length ( LL) of a single variable index model is more suited for a leaf 
area estimation model; LA = 0.1 x (1 + LL) "(R^ 2 0.77) , x? = 0:39. The results demonstrated that this is a method of 
instrument calibration and accurate estimation of other leaf traits of Chinese fir, providing data about single leaf area and 
improving the model accuracy and stability for Chinese fir leaf area. estimation. Moreover, this approach is effective in 
providing data to support advice to plantation management, and for verification and improvement of leaf morphology and leaf 


functional traits of other species. 


Key Words: Cunninghamia lanceolata; morphological characteristics ; leaf area; estimation model 


植物 叶片 被 认为 是 碳水 交换 和 能 量 平衡 过 程 最 重要 的 场所 … ,叶片 各 工作 环节 对 外 界 环境 的 响应 特点 、 
功能 发 挥 .协调 机 制 和 权 衔 策略 ,以 及 其 结构 建成 和 形态 变异 成 为 近年 来 的 研究 热点 。 目 前 对 叶片 结构 的 研 
究 已 经 从 单 株 植物 的 叶片 功能 性 状 扩展 到 群落 冠 层 、 区 域 生产 经 营 、 其 至 全 球 尺度 下 的 生物 功能 地 理 参 考量 ， 
为 全 球 生产 力 变 化 和 代谢 理论 提供 依据 和 验证 。 植 物 叶 面积 ( Leaf area,14) 是 指 植物 单 叶 或 总 平均 叶片 
的 单 面 投 影 面积 中 ;以 [下 为 基础 转换 成 的 相关 指标 ,包括 比 叶 面积 ( Specific leaf area,SLA , 单 叶 面积 / 叶 干 
重 ,cm /g)、 林 分 叶 面 积 指 数 (Leaf area. index, LAI, 叶 面 积 /土地 面积 ,mm ) ,植被 季节 性 叶 面 积 动 态 
(NDVI) 等 ; 叶 面 积 大 小 能 耳 接 影响 叶片 光 截 取 和 碳 获 取 能 力 ,是 作物 育种 农林 生产 经 营 、 模 型 佑 算 以 及 物种 
结构 变异 -功能 适应 机 理 分 析 的 关键 参考 量 。 

LA 与 植物 光合 .呼吸 蒸腾 等 生理 生态 过 程 密切 相关 ,因此 常 被 用 来 对 农作物 和 果树 进行 品种 选 育 和 性 
状 培育 ,以 及 估算 农作物 产量 和 森林 生产 及 其 影响 机 制 ,初步 探讨 生态 系统 服务 功能 ,促进 农林 生产 经 
营 。 基 于 冠 层 光合 模型 "生物 量 模 型 ”| 全球 植被 动力 学 模型 (DGVM)'"" 生物 地 球 化 学 循环 过 程 模 型 
(Biome-BGC) ”等 生态 学 模型 研究 ,1L4 与 其 相关 的 指标 都 作为 植物 重要 特征 参数 参与 模型 估算 输入 与 输 
出 。 还 有 研究 通过 对 叶 化 石 的 LA. 等 形态 学 指标 测量 来 重建 古 气候 历史 和 叶 的 进化 ,以 及 对 叶 形 态 指 标 
单独 或 联合 研究 来 解释 个 体 在 不 同 环境 下 的 可 塑性 和 适应 机 制 ……” ,以 此 来 分 析 物 种 结构 变异 -功能 适应 机 
理 075€. 但 SLA 和 LAI 的 值 受 叶 大 小 的 影响 较 大 一 ” , 叶 面 积 准 确 测 定 是 研究 植物 光合 生产 和 生理 生态 
等 过 程 的 基础 ,因此 ,精确 的 LA. 值 及 其 测量 具有 重要 意义 。 

杉木 (Cunninghamia lanceolata.) 作为 我 国 分 布 最 为 广泛 的 造林 树种 ,对 杉木 研究 主要 集中 在 与 林 分 生 
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产 和 生物 量 之 间 的 关系 “” 。 杉 木 叶片 呈 披 针 状 , 单 叶 小 而 多 。 已 有 报道 的 通过 仪器 间接 测量 的 松 类 LA 处 
于 2.242 一 4.364 cm?" , 仅 有 研究 实测 杉木 单 叶 LA 处 于 0.6873 一 1.316 cm? 7*5, 杉木 SLA 通过 模型 预测 处 
于 37 一 50 em^/g. ^ ^* ,手工 实测 值 为 79.1 em^/g 7 。 而 被 报道 的 杉木 LAI 实测 和 模型 估算 值 都 在 0.02 一 
6.34 m2]m2? 之 间 ! ”0 232 ,通过 遥感 数据 模型 估算 的 NDVI 信息 值 在 0.16 一 0.71 Z0, 42K LA 值 的 不 确 
定性 大 ,有 碍 于 生产 经 营 高 效 性 风险 评 们 有效 性 及 相关 科学 人 研究 进展 ,其 LA 的 测量 难度 大 是 造成 其 值 的 不 
确定 性 的 原因 之 一 。 

LA 测量 可 以 分 为 人 工 测 量 法 和 仪 带 测量 法 。 人 工 测 量 法 消耗 人 力 物 力 较 多 ,不 适用 于 大 批量 的 1L4 测 
定 , 但 结果 可 靠 ,和 常用 于 检验 和 校准 间接 测量 结果 。 仪 器 测量 法 对 样本 有 选择 性 ” ,日 受 叶 片 自身 平整 度 
以 及 测量 视角 的 不 同 而 影响 精度 ” , 旦 若 样 品 残 缺 较 大 ,也 会 产生 较 大 的 误差 。 目 前 国内 外 地 面 实测 工作 大 
多 采用 仪 锅 测量 方法 ,缺乏 人 工 测 量 方法 校 验 ,不 同方 法 之 间 缺 乏 比 较 。 友 的 测定 也 可 基于 简单 易 测量 的 指 
标 叶 长 和 叶 宽 ,或 是 二 者 的 综合 指标 如 长 宽 比 、 周 长 等 进行 回归 模型 预测 ” ,已 有 大 量 研究 进行 广泛 验证 ,并 
认为 该 方法 高 效 快捷 。 预 测 模型 一 般 仪 选用 单 变 量 模 型 或 多 变量 模型 ,对 两 种 模型 同时 进行 比较 人 研究 较 少 ， 
昌 样 本 大 多 集中 在 农产品 ”果树 、 花 开 生 全 等 ,以 及 阔 叶 常 绿 乔 木 ”等 ;它们 都 属于 益 叶 规则 叶 树 
种 ,相对 针叶树 种 而 言 叶 形 态 指标 测定 更 简便 。 目 前 基于 叶 形 态 指标 对 针叶树 种 建立 叶 面积 估算 模型 的 研究 
人 鲜 有 报道 。 

本 人 研究 对 杉木 LA 测定 采用 游标 卡尺 法 和 手持 叶 面 积 仪 法 对 杉木 单 叶 LA 进行 同时 测量 ,得 到 叶 长 ( Leaf 
length, LL) 平均 叶 宽 (Mean leaf width, LW )、 最 大 叶 完 (Max leaf width; LW a) 平均 叶 厚 (Mean leaf 
thickness, LT,..,) .最 大 叶 厚 ( Max leaf thickness, LT...) 、 叶 面积 ( Leaf area, LA) , 叶 周 长 (Leaf perimeter, LP) , 
叶 延 长 率 (Leaf elongation, LE) 8 个 形态 学 指标 ,分 别 通过 多 变量 线性 模型 和 单 变量 非 线 性 模型 构建 杉木 LA 
估算 模型 。 无 论 对 校准 仪器 测量 结果 还 是 利用 杉木 LA. 准确 估 测 其 他 叶 功 能 性 状 指 标 ,以 及 促进 相关 的 连续 
性 研究 都 具有 现实 意义 ,并 且 补 充 杉木 单 叶 面 积 数据 空 自 ,提高 杉木 叶 面 积 佑 算 模 型 精度 及 稳定 性 ,为 今后 杉 
木 树种 研究 及 提高 杉木 人 工 林 经 营 管理 提供 基础 数据 ,为 验证 和 完善 其 他 树种 的 叶 卢 性 状 及 其 叶 所 功能 性 状 
的 测定 提供 借鉴 意义 。 


1 材料 与 方法 


1.1 研究 地 概况 
本 研究 试验 设 在 中 国 杉木 中 心 产 区 的 湖南 会 同 杉 木林 生态 系统 国家 野外 科学 观测 研究 站 (109?"45' 工 ， 
26°50'N) 和 非 中 心 产 区 的 河南 信阳 南 湾 实 验 林场 (113°58' 玉 ,31°53'N) ,其 海拔 高 度 范 围 分 别 为 270 一 400、 
300—1100 m, 均 系 低 山地 丘陵 地 貌 。 湖 南 会 同 与 河南 信阳 同属 于 亚热带 季风 气候 ,气候 温和 ,年 均 温 在 
15.1 一 16.8% 间 ; 年 均 相 对 湿度 均 在 77% 以 上 ,十 热 同 期 ,年 均 降 水 量 分 别 为 1304.2 mm 和 1106 mm, 雨水 集中 
在 夏季 ,日照 充沛 ,年 日 照 时 数 均 在 1900 h 以 上 。 湖 南 会 同 林 分 密度 约 为 2400 株 /hm- ,河南 信阳 约 为 1500— 
1800 P&/hm^ ,上 且 士 壤 类 型 均 为 山地 黄 壤 , 其 光照 和 水 热 条 件 均 适 宜 杉木 生长 。 
1.2. .sS]2E 2r i E 
于 杉木 生长 旺盛 期 4 月 在 湖南 会 同和 河南 信阳 的 T 
不 同 林 龄 ( 幼 龄 林 、 中 龄 林 、 近 熟 林 、 成 熟 林 ) 中 设置 样 
Jr (10-mx10 m) ,LA 主要 受 林 分 的 光照 和 水 分 等 多 重 
影响 “… ,因此 在 每 个 样 方 内 随机 选取 3 株 杉 木 ,为 保 
持 叶 龄 序列 ,分 别 在 树冠 上 层 中层、 下 层 随机 采集 生长 图 1 用 于 形态 指标 测定 的 杉木 小 时 示例 
健 康 . X 病 "n 害 的 $ j Ad i B 为 防 En 十 片 失 水 | 将 采 集 we of leaf morphological traits measurements of 
样品 尽快 带 回 实验 室 。 由 于 杉木 叶 形 不 是 常规 的 针 叶 
JE. ,将 每 束 针 叶 近 似 假 设 成 圆柱 体 对 针叶树 种 进行 叶 面 积 直接 测量 的 方法 并 不 适用 一 。 我 们 选择 对 杉木 
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每 片 单 叶 进 行 逐 个 测量 ,在 每 个 枝条 上 随机 选择 三 片 小 叶 ,逐一 编号 并 通过 两 种 方法 分 别 测 定 小 叶 形 态 指标 
(图 1)。(1) 游 标 卡 尺 法 ;用 游标 卡尺 (精确 到 0.002 cem) 分别 测 量 杉木 小 叶 的 LL( 叶 连接 叶柄 的 最 底 端 到 叶 
梢 的 距离 ) ,小 叶 1/4、1/2、3/4 长 处 的 LW 和 LT” ,每 处 重复 3 次 取 平 均值 作为 该 处 的 LW 和 ZL7, 计 算 3 处 平 
均值 的 均值 作为 该 叶片 的 LW 和 ZT, 求 得 LL LW nean x LW max LI mean LT, Se 5 个 直接 测量 指标 。 将 杉木 小 叶 
的 投影 面积 近似 看 做 和 矩 形 ,将 LL 与 LW 的 乘积 算得 LA ,LL 与 LW 和 的 两 倍 为 LP,LL/LW, 算 得 LE, 得 到 
3 个 间接 转 算 指标 。(2) 手 持 叶 面积 仪 法 (Yaxin- 1242) ,得 到 LL( 精确 到 0. 1em) LW nean MR] 0.01 em) 56 2 
个 直接 测量 指标 ,以 及 LA( 精 确 到 0.001 em?) .LP( 精 确 到 0.01 em) .LE( 保 留 三 位 小 数 )3 个 间接 转 算 指 标 。 
1.3 ”模型 选择 和 拟 合 

为 了 构建 不 同形 态 指标 与 LA 最 佳 拟 合 估算 模型 ,采用 多 变量 线性 回归 模型 和 非 线 性 回归 指数 模型 进行 
模型 拟 合 。 首 先 , 为 了 探索 不 同 指标 对 LA 的 影响 和 贡献 量 , 选 择 综合 指标 共同 作用 的 多 变量 线性 模型 ;Y=at+t 
Èb X, i=1 一 7, 并 且 通 过 逐步 法 对 变量 第 选 进行 模型 优化 。 为 了 探究 不 同 指标 单独 对 LA 的 估算 ,根据 散 点 
图 分 布 情况 ,采用 单 变量 非 线 性 回归 指数 模型 .LA4= a x(1 + X)" ,通过 对 校正 的 决定 系数 (R )\ 卡 方 (x ) 和 
显著 性 (PP) 综合 参考 ,对 变量 拟 合 和 研究 ,并 且 筛 选 出 最 佳 拟 合 指标 。 
1.4 数据 分 析 

数据 进行 初步 的 整理 在 Microsoft Excel 2010 上 进行 ,数据 分 析 和 模型 拟 合 用 IBM SPSS Statistics 21.0 , 作 
图 在 R 软件 上 完成 。 


2 结果 与 分 析 


2.1 杉木 叶 形 态 指 标的 生态 特征 
Æ LL LW, LW, LT LT, .LALE 和 LP 等 8 


mean 80 
个 指标 实测 值 的 分 布 中 ( 表 1) ,变异 系数 大 小 为 14> 
LE SLT un >LE>LL>LP>LW mean >LW x ,其 中 LA 的 灾民 D 60 
fæ 
程度 最 大 ( CV=0.513) ,均值 为 (0.797+0.409) (SD)em? È 
(图 2) ,数值 大 致 分 布 在 (0.758 一 0.836) cm2(95% CI) 。 g 40 
R 


LT ax LT, PA TR ER E Se BE DUX-T. LA(CV 5 0.489, 
0.465) ,LT,.. I LT... 均值 分 别 为 (0.032+0.016) cm 和 20 
(0.028+0.013) cm, 但 数值 分 布 集中 呈现 明 显 的 尖峰 态 


(IE HE 4) 8I 16.383 81 13.851) ,数据 分 布 与 均值 最 接 


0 
0 02 04 06 08 10 12 14 16 18 20 22 


Vr CLT,.. 8 LT ag MAD.SE ,9596CI 分 别 0.011 0.001 , 叶 面 积 Leaf area/cm? 
0.03 一 0.034 和 0.009.,0.001,0.027—0.029) , LE, LL, H2 EGERSSRNAG 
LP3 个 指标 变异 系数 接近 ( CV=0.414, 0.411, 0.398) , Fig.2 Frequency distribution of leaf area of C. lanceolata 


Hp LE 数值 分 布 范 围 最 大 且 值 偏差 最 大 (0.526— 
33.495, MAD=4.628 , SE 2 0.292 ,9596 CI y 13.941—15.087) , LL 和 LP 的 数值 分 布 总 体 分 布 稳定 ,和 平均 值 相 
差 不 大 ( MAD 分 别 为 1.071 和 2.19) 。 变 异 程度 最 小 为 衡量 叶 宽 的 两 个 指标 LW,. Fl LW aa ( CV 2 0.238, 
0.224) ,二 者 数值 波动 和 偏差 都 较 小 ,分 布 集中 ,1 到 和 LW, 的 范围 .MAD SE 和 9596 CI 分别 为 0.11 一 
0.486 cm 和 0.8 一 0.377 cm ,0.051 和 0.043 ,0.003 和 0.003 ,0.284 一 0.297 cm 和 0.227 一 0.238 cm。 
2.0 形态 学 指标 相关 性 

通过 对 杉木 叶 形 态 学 指标 相关 性 (Pearson) 分 析 ( 表 2) 发 现 不 同 指标 对 LA 的 直接 影响 或 间接 影响 。Z4 
除了 与 LE LP 线性 相关 程度 较 弱 (r= 0.292,0.244) ,与 其 他 指标 都 呈 不 同 程度 的 显著 线性 正 相 关 性 。 其 中 LA 
与 LL LT... ERRE EIA% (720.896, 0.911) ,LL 5 LT ， 也 成 极 显著 线性 正 相 关 (r=0.999) ; LA 与 LW,.,、 
LT ， 呈 较 强 线性 正 相 关 (r= 0.682,， 0.709) ,LW 与 L7, 之 间 也 呈 极 强 线性 正 相 关 (r=0.94)。 因 此 ,1L4 与 


http ://www.ecologica.cn 


10 期 Wk  SEIRLZKRUDJEGSREAE S UT TRURRA EC 3573 


LT pean LT BEBE TH F LL II LW pean IRAR BRI LW oa 5 LE LP 几乎 不 呈 相 关 性 (7 分别 为 0.073 
和 0.041) ,其 余 指标 间 都 呈 不 同 程度 的 显著 相关 性 。LW,, Ej LW nean LT pa ATAR DR] 5 BUE — F6 f 4 2: , fH. 
相关 程度 较 弱 (r=-0.29, r2 0.201) , LE LP 与 其 他 指标 之 间 线 性 相关 都 较 弱 (r<3.5) ,但 二 者 之 间 呈 极 显 
著 的 线性 正 相 关 (r=0.946 ) 。 


R1 杉木 叶片 形态 指标 值 总 体 分 布 特征 


Table 1 Distribution characteristics of C. lanceolata leaf morphometric value 


统计 量 指标 Measurements 
Statistics LA/ cm? LL/cm LW cm LW nax €M LT aean” cni LT, ,,/ cm LE LP/cm 
样本 量 N total 424 424 424 385 385 385 424 424 
均值 Mean 0.797 3.344 0.233 0.29 0.028 0.032 14:514 7.097 
标准 差 n 0.409 1.375 0.055 0.065 0.013 0.016 6.003 2.825 
Standard Deviation 
J^ L X] Y». c p 
D MEETS 0.02 0.067 0.003 0.003 0.001 0.001 0.292 0.137 
SE of Mean 
均值 的 95% 置 信 区 间 上 限 
Lower 9595CI of Mean 0.758 3.213 0.227 0.284 0.027 0.03 13.941 6.8277 
均值 的 95% 置 信 区 间 下 限 

0.836 3.475 0.238 0.297 0.029 0.034 15.087 7.366 
Upper 9596 CI of Mean 
偏 度 Skewness 0.214 -0.02 -0.33 -0.104 2.143 2.721 0.443 -0.12 
峰 度 Kurtosis —0.26 0.099 0.064 0.099 13.851 16.38 0.459 0.138 
dE ark 
变异 系数 MEV 0.513 0.411 0.238 0.224 0.465 0.489 0.414 0.398 
Coefficient of Variation 
SZ T^ a s 
FIENT 2E TN 0.327 1.071 0.043 0.051 0.009 0.011 4.628 2.19 
Mean absolute Deviation 
最 小 值 Minimum 0.017 0.1 0.08 0.11 0.003 0.004 0.526 0.17 
最 大 值 Maximum 2.166 8.332 0.377 0.486 0.137 0.14 33.495 17.184 


iE: 14: 叶 面积 ,leaf area; LL: 叶 长 ,leaf length; LW :平均 叶 宽 ,leaf mean width; LW pas: E HT 95, leaf maximum width; LT pean :平均 叶 厚 ， 


mean max 


leaf mean thickness; LT,.. :最 大 叶 厚 ,leaf maximum width; LE : IHE K% leaf elongation; LP : HT J&]& , leaf perimeter 


表 2 杉木 叶片 形态 学 指标 相关 性 (Pearson ) 分 析 


Table 2 One-way analysis of variance (Pearson) of Chinese fir leaf morphometric 


指标 Measurements LA/cm? LL/cm LW, [cm LW,,/em LT /em LT /em LE TBI 
LA/cm? 

LL/cm 0.896 ** 

LW /em 0:682." xus ne 

LW /em gas 0.787 ** ET 

LT, / cm 0.9117* 0.999 ** 0.359 ** 0.764 ** 

LT, em 0.709 ** 0.384 ** 0.905** — -0.201** 0.417** 

LE 0.292 ** ]4 0.288 ** 0.073 0.248 ** 0.344 ** 

LP/cm 0.244 ** 0.19** 0.253 ** 0.041 0.197 ** ge 0.946 ** 


* x 在 P=0:05 水 平 上 表现 显著 性 差异 。LA: 叶 面积 ,leaf area; LL; WE T , leaf length; LW pean :平均 叶 宽 ,leaf mean width; 工矿， :最 大 叶 宽 ， 
leaf maximum width; LT ean: F H5 HF FÉ, leaf mean thickness; LT,,. : EE X HT JE, leaf maximum width; LE: 叶 延长 率 , leaf elongation; LP; 叶 周 长 ， 


leaf perimeter 


2.3 杉木 叶 面 积 多 元 回归 模型 

将 LL LW a LW n LT a LT, LE LP 等 7 个 指标 与 LA 进行 多 变量 线性 拟 合 和 优化 ,由 上 述 杉木 叶 
形态 学 指标 相关 性 (Pearson ) 分 析 ( 表 2) 得 出 LL 5j LT, 存在 的 极 强 显著 相关 性 (r=0.999) , HKF LL 5 LA 
的 相关 性 (r=0.896) ,这 违背 了 变量 间 的 相互 独立 原则 ,为 了 避免 多 重 共 线性 市 来 的 影响 ,在 多 变量 线性 回归 
时 ,指标 LL 被 排除 。 
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如 表 3 ,首先 ,选择 全 和 法 构建 多 变量 线性 回归 方程 模型 1: 

Y- -0.397+0.164X -0.022X, -1.253X,*0.24X,-1.416X,*0.854X,( R^ 20.981, SE =0.053) ,其 中 XX 一 XX 分 
3I LP LE LW gas W aas LT onas LT mo 

虽然 模型 拟 合 度 很 好 ( R* 20.981) ,但 存在 3 个 回归 系数 不 显著 的 指标 LW, L aean a Laan CP 70.074, 
0.031, 0.11) 。 为 使 方程 达到 最 优 ,选择 逐步 法 再 次 构建 多 变量 线性 回归 方程 。 模 型 2 以 LP 为 自 变 量 得 到 | 
Y=-0.174+0.138X (R^ 20.83, SEZ0.157) ;继续 引入 变量 LE 得 到 模型 3:7= -0.076+0.207X -0.041X, (R^ — 
0.977, SE=0.058) , X, X IIIN LP 和 LE; 进一步 引入 变量 LW 得 到 模型 4.Y= 一 0.388+0.165X| -0:023X,+ 
1.453X,( R° 20.981, SE-0.053) ,X,—X,2Y 929 LP LE LW ,此 时 LW LT, LT, 由 于 指标 不 显著 被 排 


除 (P>0.01) ,模型 4 以 最 少 的 变量 达到 最 大 的 拟 合 优 度 (R =0.981) ,模型 达到 最 优 。 


表 3 杉木 形态 学 指标 多 变量 线性 回归 
Table 3 Multiple linear regression of Chinese fir leaf morphometric 


HX IB 
系数 Coefficient 95% 置信 区 间 


方法 yia > EOS pe ER xt 95% of Confidence Interval 
Methodi G un X 系数 全 au Se TUE 
Value p Lower Upper 
ZA 模型 1 (常量 ) -0.397 0.035 0 -0.467 -0.327 
Complete subsitituion R? 2 0.981 , 标准 LP( cm) 0.164 0.005 0 0.154 0.173 
method 估计 的 误差 =0.053 — LE -0.022 0.002. 0 -0.026 -0.018 
LW mean ecm) 1.253 0.213 0 0.834 1.671 
LW max( €M) 0.24 0.134 0.074 -0.023 0.504 
LT nemn €M) -1.416 0.654 0.031 -2.701 -0.131 
LT nax( cm) 0:854 0.533 0.11 -0.195 1.902 
逐步 模型 2 (常量 ) -0.174 0.025 0 -0.223 -0.125 
RIEN OD pue Lf (e) 0.138 0.003 0 0.132 0.144 
模型 3 (常量 ) -0.076 0.009 0 -0.094 -0.057 
R? =0.977, 标准 LP ( cm) 0.207 0.002 0 0.203 0.21 
估计 的 误差 =0.058” LE -0.041 0.001 0 -0.042 -0.039 
模型 4 (常量 ) -0.388 0.036 0 -0.458 -0.318 
R? =0:981, 标准 LP( cm) 0.165 0.005 0 0.156 0.175 
估计 的 误差 =0.053 LE —0.023 0.002 0 -0.027 -0.018 
LW aean (cm) 1.453 0.16 0 1.138 1.769 


2.4 ”杉木 时 面积 指数 模型 

通过 多 元 线性 模型 发 现 不 同 指标 对 LA. 的 综合 影响 和 不 同 页 献 量 后 ,为 了 进一步 探究 单 变量 对 LA 的 信 
算 ; 根 据 不 同 指标 散 点 图 分 布 情况 ,排除 线性 模型 ,而 采用 非 线 性 模型 [4= a (1 + X) 对 7 个 形态 学 指标 和 
LA 分 别 进行 模型 拟 合 ( 表 4) 。 所 有 指标 与 LA 都 表现 出 极 强 显著 性 (P=0) ,7 个 形态 学 指标 都 对 LA 有 不 同 
程度 的 关联 和 影响 ,可 以 建立 模型 表达 它们 之 间 的 关系 。R 从 大 到 小 依次 为 :LP>LL>LW,,,>LW,,>LE>LT,, 


SLT X. 从 小 到 大 依次 为 :LP<LL<LW, <LW, SLE <ET poaa ELT, , RU 大 小 所 体现 出 的 拟 合 优 度 基本 一 


Sic. 因此 综合 二 者 比较 , 拟 合 优 度 最 好 的 是 LP.LA=0.046x(1+LP)'(R*=0.825,X 20.029) , 且 在 95% 似 然 区 
间 (95%PL) 拟 合 效果 最 好 (图 3) , RÆ LL:LAS0.1x (1-LL) 95( R*=0.77,x?=0.39) ,其 模型 拟 合 度 与 95% 
PL 的 拟 合 度 都 与 LP 相似 。 而 LW, E; LW, 拟 合 效果 一 般 , 模 型 分 别 为 :LA=0.144x(1+LW,.,,)” (R^ = 


0.504,x 20.083) , LA- 0.18x ( 1L -LW,.. 9 ?9( R? 20.484,47 =0.075)。L7T, 和 LT 两 个 模型 拟 合 效果 较 差 


max mean max 


( R? =0.062,0.042) ,不 建议 作为 LA 估算 模型 。 
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表 4 杉木 叶 面 积 指数 模型 


Table 4 Exponential model of Chinese fir leaf morphometric 


白 恋 量 区 参数 Parameters Statistics 统计 量 

a b N R? y P 
叶 长 (LL) /em 0.1 1.398 424 0.77 0.039 0 
平均 叶 宽 (LW ) /em 0.144 7.973 424 0.504 0.083 0 
EKRI SE (LW aa) /em 0.18 5.963 385 0.484 0.075 0 
平均 叶 厚 (LT hean) /em 0.702. 6.915 385 0.062 0.136 0 
最 大 叶 厚 (LT...) /em 0.731 4.833 385 0.042 0.139 0 
叶 延 长 率 (LE) 0.153 0.61 424 0.281 0.12 0 
叶 周 长 (LP) /em 0.046 1.35 424 0.825 0.029 0 


N= Number of Points; R? = Adj; R-Square; y? =Reduced Chi-Sqr 


2.5 杉木 叶 面 积 估算 模型 检验 

通过 对 最 优 多 元 线性 模型 4: 了 = -0.388+0.165X, 一 0.023X,+1.453X,(R?* 20.981, SE-0.053, 了 为 叶 面 积 ， 
一; 分别 为 LPLE LW) ÆFA K LAS 0.046x ( 1€LP) 5( 82 20.825,43? =.0.029), 和 叶 长 :LA4=0.1x(1+ 
LL) "(R^ 20.77 x. =0.39) 的 叶 面 积 非 线性 模型 的 检验 。 多 元 线性 模型 拟 合 优 度 较 高 ,通过 LA 实测 值 和 预 
测 值 直接 的 关系 可 以 发 现 拟 合 线 接 近 Y=XX 的 拟 合 线 ( 图 4) ,实测 值 和 预测 值 接近 程度 很 高 。 由 图 4 可 以 进 
一 步 发 现 , 当 用 多 元 线性 模型 对 叶 面 积 进 行 估算 时 ,估算 误差 可 能 产生 在 叶片 极 小 (LA<0.4 cm?) 和 极 大 (7L4> 
1.6 cm ) 处 ,此 时 预测 值 小 于 实测 值 。 而 当 用 单 变量 非 线性 模型 估算 时 ,无 论 是 基于 LP s LL, —d ER EM 
算 残 差 图 皆 呈 随机 分 布 (图 5) ,分 布 趋势 接近 ,估算 误差 可 能 较 多 产生 于 LP RII LL EK (LP»12 cm 和 LL>6 
cm) 时 ,此 时 残 差 较 大 , 且 较 多 分 布 在 残 差 负 值 域 ,预测 值 俩 大 。3 个 模型 通过 检验 符合 估算 要 求 ,都 达到 恨 好 
的 拟 合 效 果 , 仅 在 叶片 极 小 或 极 大 的 情况 时 出 现 些许 误差 ,由 2.1 分 析 认 为 ,这 两 种 叶片 极 小 和 极 大 的 情况 丝 
超过 了 正常 杉木 叶 面积 范围 。 


3 讨论 


3.1 杉木 叶 面 积 特征 

相对 其 他 形态 学 指标 ,杉木 的 LA 变异 (CV=0.513) 最 大 ,这 与 对 南方 四 种 阔叶树 种 的 形态 学 指标 变异 分 
析 的 结果 相同 '” ,原因 可 能 由 于 4 度量 单位 为 平方 厘米 ,因此 测量 方法 导致 的 误差 相对 其 他 指标 贡献 更 
大 。 且 杉木 单 叶 为 狭长 的 披 针 状 , 其 LL 的 变异 程度 较 LW 更 大 (CV=0.411, 0.238) ,通过 相关 性 (Pearson ) 分 
析 ( 表 2) 发 现 LL 也 为 L4 的 变异 提供 了 更 多 的 解释 空间 (r= 0.896) ,这 与 已 有 对 阔 叶 的 人 研究 发 现 不 同 ,已 有 
研究 认为 阔 叶 的 LWX LA 的 变异 程度 有 近 90% 的 解释 力度 , 比 LL(75%) 更 大 。 可 见 不 同 叶 形 下 ,LL 和 
LW 对 L4 的 影响 和 贡献 量 不 同 。 已 有 研究 对 杉木 的 LA. 的 测量 范围 在 0.6873 一 1.316 em ,标准 差分 布 在 
0.136—0.316 cm^ zz [H] 7^! ,本 人 研究 杉木 LA 的 95%CI 范围 在 0.758 一 0.836 ecm’, 而 本 研究 的 有 效 值 的 范围 相 
比 已 有 研究 更 精确 ,但 标准 偏差 偏 大 (SD=0.409) ,原因 可 能 是 由 于 已 有 研究 的 样本 主要 是 一 棵 树 的 50 一 60 
杉林 叶片 ,而 本 研究 的 样本 在 不 同 林 龄 .不同 冠 层 和 不 同 叶 龄 丝 有 分 布 , 跨 度 相 对 已 有 人 研究 更 广 , 且 和 手工 测量 
误差 较 小 ,因此 有 效 值 更 精确 ,而 标准 偏差 较 大 。LT 较 大 的 变异 度 与 集中 的 数值 分 布 体现 了 杉木 叶片 L7 DU 
量 的 困难 ,要 求 精确 度 高 ,0.01 em 带 来 的 影响 相 比 其 他 指标 更 明显 。 但 由 于 杉木 叶片 厚度 测量 只 由 人 工 游 标 
卡尺 测量 ,精确 度 得 到 保障 。LE 的 数值 偏差 大 是 由 于 杉木 叶 较 长 而 狭窄 ,而 LE 为 通过 测量 LL 和 LW aean FEA 
而 来 的 比 指标 ,因此 由 于 测量 产生 的 误差 会 比较 显著 ,导致 LE 波动 范围 更 大 。LL、LP LW 的 数值 总 体 分 布 集 
中 ,偏差 小 ,可 能 原因 为 杉木 单 叶 较 小 ,自身 LL 和 LW 相对 其 他 树种 变异 小 , 且 测 量 方法 精确 ,测量 方法 带 来 
的 总 体 误差 小 ,这 为 模型 的 建立 提供 依据 。 

LE 与 LP 为 非 直 接 测 量 的 转 算 指标 ,而 LE 和 LP 的 转 算 方式 与 LW 没有 直接 关系 ,因此 LW, 与 LE LP 
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0.5 
Ü 
0 1 2 3 4 5 6 7 8 9 0 0.1 0.2 0.3 0.4 0.5 
mHE LL 平均 叶 宽 LW mean 
Leaf length/cm Mean leaf width/cm 


0 0.1 0.2 0.3 0.4 0.5 0 0.05 0.10 0.15 
RARI RE LW max 平均 叶 厚 LT mean 
Max. leaf width/cm Mean leaf thickness/cm 


叶 面 积 LA 
Leaf area/ cm? 


0 0.05 0.10 0.15 
最 大 叶 厚 LT max 


Max. leaf thickness/cm Leaf elongation 


— FitofLA 
— 95% CL 
95% PL 


0 5 10 15 20 
叶 周 长 LP 


Leaf perimeter/cm 


图 3 杉木 叶 面积 指数 非 线 性 拟 合 
Fig.3 Nonlinear fitting of Chinese fir leaf area 


Fit of LA : 叶 面 积 拟 合 线 ,Fit of leaf area;95% CL:95% 置 信 区 间 线 ,95% confidential line;95% PL:95% probability line 


几乎 不 呈 相 关 性 。 很 明显 ,Lh 会 随 着 LL 或 LW 的 增 大 而 增 大 ,这 与 对 常 绿 阔叶树 种 山 白兰 的 研究 相同 , 
但 LT 、LT,, 可 能 是 由 于 LL 和 LW 的 影响 而 呈 间 接 相 关 。LE 和 LP 都 是 用 相同 的 两 个 指标 (LL、LW,.,) 
转 算 而 来 ,只 是 通过 不 同 的 运算 方式 ,因此 二 者 有 明显 的 线性 相关 性 。 
3.2 ”杉木 叶 面 积 估算 模型 

由 多 变量 线性 回归 方程 发 现 LPLE LW 对 74 的 直接 影响 最 大 上 是 最 显著 ,而 由 于 ZL 在 模型 建立 前 被 排 
除 , 但 LP 和 LE 都 为 LL 和 LW 的 间接 转 算 指 标 ,因此 可 以 认为 LL 和 LW 都 对 LA 的 直接 显著 影响 ,已 知 LL 
fl LW 的 值 , 即 可 以 通过 转 算出 LE 和 LP 对 叶 面 积 进 行 估算 。 存 已 有 研究 中 ,LL 和 LW 也 经 常 作为 综合 指标 
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共同 对 叶 面 积 进行 估算 5 ” ”1 。1.9% 未 拟 合 部 分 " 
造成 的 可 能 原因 是 操作 方法 带 来 的 误差 ,是 可 以 接 
5E RJ s ij PA evt 
<.5 14 fil 

多 变量 线性 模型 可 以 建立 高 精确 度 的 LA 预测 模 EE 12 
型 ,但 其 需要 多 个 变量 计算 往往 会 造成 计算 的 复杂 和 宛 Bud 
余 , 从 简便 性 上 考虑 ,本 人 研究 进一步 建立 单 变 量 非 线性 du n. . 5 4 Y= -0.01984 + 1.01201X 
模型 对 LA 进行 拟 合 。 在 非 线性 拟 合 模型 中 ,IP RIL 3 02 L at Wi 

` m ^ HM 0 do M 

的 拟 合 较 好 ( R* 20.825, 422 0.029) ,这 与 已 有 研究 对 -0.2 BS > 
E ; T LX Fio .a [54-56] EART : 
ISI TERURLIO a REARUR 5x b EE HK 0 02 04 06 08 10 12 14 16 1.8 20 22 
际 考虑 ， 由 于 LP 在 测量 中 的 相对 不 简便 性 1 ,更 建议 叶 面 积 Leaf area/cm? 


LL 来 对 杉木 小 叶 LA 进行 模型 :LA=0.1x(1+ 

a (It qa 杉木 叶 面积 实测 信 与 表 3 中 多 元 线性 身 内 模型 4 的 估算 人 
LL)" (R 20.77, 六 =0.039) , 且 在 拟 合 中 的 部 分 剩余 yez 

的 2396 AU 2s 可 能 由 于 测量 带 来 的 误差 , DE Fig.4 Relationship between measured and estimated leaf area 
E AHE E i Ai [E709 , 由 此 可 见 测 量 方法 的 精确 度 在 when leaf area was modeled implicitly with leaf perimeter, leaf 
对 LA 的 估算 中 影 响 不 J 小 视 TE 其 他 [js] IH- 物种 的 叶 面 elongation and leaf mean width ( multiple linear regression model 
积 模型 估算 中 ,LW A5] THO iced ot seem ,但 “wee? 

在 杉木 叶 面 积 研 究 中 ,LW 模型 拟 合 度 表 现 没 有 LL 好 ， 
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图 5 基于 时 周 长 和 叶 长 的 叶 面积 非 线性 拟 合 模型 残 差 图 


Fig.5 Residual of leaf area based on leaf perimeter and leaf length nonlinear fitting models 


这 与 指标 间 的 相关 性 (Pearson ) 219r (K 2) 得 出 的 结论 相同 。 是 因为 杉木 叶 形 态 窄 而 狭长 ,以 em 为 度量 单 
位 ,LW 数据 样本 分 布 集中 ,0.01 em 的 测量 误差 带 来 的 影响 对 LW 更 大 。LE 作为 间接 转 算 指标 ,也 存在 实际 
测量 的 不 简便 性 。LT 数据 分 布 偏向 拟 合 网 的 一 边 ,原因 认为 是 叶 厚 的 所 用 度量 单位 是 厘米 所 致 ,用 毫米 可 能 
会 更 合适 ,日 95% 署 信 区 间 在 叶 厚 大 于 0.05 em 处 变 大 ,原因 推测 认为 杉木 叶 厚 普遍 小 于 0.05 em, KF 0.05 
cm 的 样本 不 多 , 且 个 别离 奇 值 的 出 现 可 能 原因 推测 为 由 于 测量 叶 厚 的 难度 大 且 要 求 精度 高 ,手工 测量 造成 的 
误差 所 致 ,由 与 LT LT 可 能 是 由 于 友和 7 的 影响 而 成 伪 相 关 , 这 与 现实 中 LA 5 LT 为 间接 影响 观 
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SBU ,其 直接 影响 不 明显 ,需要 更 进一步 的 研究 。 


4 结论 


本 研究 通过 对 LA 与 7 个 相关 的 叶 形 态 测量 指标 (LLW,, LW max LT nean LT LE ZLP) 间 的 相关 分 析 和 
模型 拟 合 ,杉木 叶 面积 为 (0.797+0.409) cm ,大 部 分 分 布 在 0.758—0.836 cm ,不 同 的 指标 对 LA 的 直接 或 间 
gzsé pA [8] , EP LL AI LW xt LA 的 直接 影响 更 大 ,LT 对 LA 以 间接 影响 为 主 。 通 过 进一步 基于 叶 形 态 指标 对 
LA 进行 拟 合 估算 ,从 简便 性 上 考虑 ,LL 的 单 变量 指数 模型 适合 对 LA 进行 估算 LA4=0.1x(1+LL)(R*=0.77,， 
x =0.39) 。 从 精确 角度 考虑 ,LL 和 LW 的 多 元 线性 模型 更 好 ,Y= 一 0.388+0.165X, 一 0.023X,+1.453X,( R= 
0.981,8E 20.053) ,Xi 一 X; 分 别 为 LP、LE LW。 通过 模型 检验 后 认为 两 种 模型 丝 符 合 叶 面积 估算 要 求 ; 飘 
型 较 大 误差 情况 可 能 出 现在 叶片 极 小 或 极 大 时 ,但 该 值 域 已 超过 正常 杉木 叶 面 积 范 围 。 根 据 不 同 的 情况 选择 
不 同 的 模型 ,能 对 杉木 叶 面 积 进行 非 破坏 性 测量 , 既 避 免 了 仪器 带 来 的 限制 ,也 更 加 省 时 省 力 ; 为 今后 杉木 的 
叶 功 能 性 状 的 测量 提供 更 大 的 简便 性 ,为 今后 相关 的 连续 性 研究 提供 参考 和 借鉴 意义 - 
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